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Abstract 

Solutions to the equation dt£(x,t) — ^A£(x,t) = X\S(x, t)\ 2 £(x, t) are investigated, where 
S(x, t) is a complex Gaussian field with zero mean and specified covariance, and m ^ is a complex 
mass with Im(m) > 0. For real m this equation describes the backscattering of a smoothed laser 
beam by an optically active medium. Assuming that S(x, t) is the sum of a finite number of 
independent complex Gaussian random variables, we obtain an expression for the value of A at 
which the g-th moment of \£ (x,t)\ w.r.t. the Gaussian field S diverges. This value is found to be 
less or equal for all Im(m) > and \m\ < +oo than for |m| = +oo, i.e. when the A£ term 

is absent. Our solution is based on a distributional formulation of the Feynman path-integral and 
the Paley- Wiener theorem. 
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I. INTRODUCTION 



We investigate the breakdown of linear amplification in a system driven by the square 
of a Gaussian noise. This problem which models the backscattering of an incoherent laser 
by an optically active medium was first considered by Akhmanov et al. in nonlinear optics 
[1], and by Rose and DuBois in laser-plasma interaction [10]. The latter investigated the 
divergence of the average solution to the stochastic PDE 

d t £(x,t) - ±A£(x,t) = X\S(x,t)\ 2 £(x,t), 
t > 0, xeAc M. d , and £(x, 0) = 1, 

heuristically and numerically in the "diffractive case" where Im(m) = and Re(m) ^ 0. 
Here A > is the coupling constant and S is a complex Gaussian noise with zero mean 1 . 
More recently, this problem was analyzed from a more rigorous mathematical point of view 
in [2] and [7]. The "diffusive case" in which Re(m) = and Im(m) > was considered in 
[2], and the one dimensional diffractive case was considered in [7] for a restrictive class of 
S"s. In the present work we will consider the general case m/0 and Im(m) > for A a 
rf-dimensional torus with d < 3. As in [2] and [7] we will express the solution to (1) formally 
as the Feynman-Kac path-integral 

£( X ,t) = [ e fo[T*W a +^(*(r),r)|"]«»r d[a . ( . )]j (2) 

Jx(-)eB(x,t) 

where B(x,t) denotes the set of all the continuous paths in A satisfying x(t) = x. 

In the diffusive case, the right-hand side of (2) is just the Wiener integral of 
exp ^A J* * \S(x(t), t)| 2 drj over B(x,t). This was used in [2] to prove, under some rea- 
sonable assumptions on the covariance of S, that for every t > and any positive integer 
q the average of £(x,t) q over the realizations of S, (£(x,t) q ), diverges as A increases past 
some critical value smaller (or equal) than in the diffusion-free case (i.e. when \m\ = +oo), 
with equality holding for a class of S. It was conjectured there that this inequality should 
also apply when diffusion is replaced by diffraction, i.e. m real, m^fl, the case of physical 
interest considered by Rose and DuBois in [10]. 

1 This is the case of interest in laser-plasma interaction and nonlinear optics in which S is the (complex) 
time-enveloppe of the laser electric field. With the help of some minor modifications, our results carry 
over straightforwardly to the cases where S is real. 
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The diffractive case is much more difficult because the right-hand side of (2) is no longer 
well defined and one cannot a priori exclude the possibility that destructive interference 
between paths makes the sum of divergent contributions finite, raising (possibly to infinity) 
the critical value of A at which the average of (2) diverges. Using heuristics and numerical 
simulations, Rose and DuBois argued that (\£(x,t)\ 2 ) should diverge for every t > as 
A increases to some finite critical value [10]. The conjecture made in [2] that diffraction 
should actually lower the critical coupling (or, at least, not increase it) compared to the case 
\m\ = +oo was proved in [7], for very special choices of S, for the divergence of (\£(x,t)\). 

In this paper, we extend the results of [7] to a much wider class of S. We analyze 
the divergence of {\S(x, t)\ q ) for any positive integer q, and we treat both the diffusive and 
diffractive cases as well as all the intermediate cases between these two limits [i.e. complex 
m with Im(ra) > and m ^ 0]. Our strategy for controlling the complex Feynman path- 
integral (2) and determining the critical value of A uses the following three ingredients: 

(i) we consider a restricted but quite wide class of S for which £(x,t) can be written as 
a Fourier-Laplace integral w.r.t. a distribution with compact support; 2 

(ii) we apply the Paley- Wiener theorem to this Fourier-Laplace integral. This yields the 
control of (2) for "large" \S\ 2 ; 

(iii) we average \£(x, t)\ q over the realizations of S and use the control obtained in (ii) to 
determine the smallest value of A for which this average blows up. 

The rest of the paper follows this strategy quite faithfully. In Section II we specify the 
class of S which we can treat. The distributional formulation of £(x,t) is given in Section 
III and the way to control its growth is explained in Section IV. Finally, the determination 
of the critical value of A and the proof of the conjecture made in [2] are given in Section V. 
It is worth noting that (i) and (ii) do not depend on S being Gaussian and thus apply also 
in a more general setting. 

2 The "Fourier-Laplace" integral w.r.t. a distribution with compact support on M. N is the continuation of 
the usual Fourier integral from to C N . 
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II. MODEL AND DEFINITIONS 



We consider the solution to the linear amplifier equation (1), written in its integral 
representation (2), with m in C+\{0}, where C+ = {m G C : Im(m) > 0}. We assume that 
S can be expressed as a finite combination of M complex Gaussian r.v., s n , 

M 

S(x,t)=J2 S n®n(x,t), (3) 
n=l 

with 



( S nS m ) — 5 n 



(4) 



The $„ are normalized such that 

(x,t)\ 2 dTd d x = 1. 



' ?1 



Furthermore, the $ n (-,r) are assumed to have second derivatives bounded uniformly in 
r G [0,t], and the <3> n (x, •) are piecewise continuous for every x G A with a finite number of 
discontinuities in [0,t] for all finite t. Note that locally, i.e. for each x and r, \S(x, r)| 2 is a 
quadratic form of 2M real Gaussian r.v. (the real and imaginary parts of the s n ), so it is a 
X 2 r.v. with 2M degrees of freedom. 

Equation (3) generalizes models of spatially smoothed laser beams in which the laser 
light is represented by a superposition of a finite number of monochromatic beamlets the 
amplitudes of which are independent r.v. [9]. For a large number of beamlets these r.v. 
can be taken as Gaussian and the laser electric field takes on the form (3) with $ n (x,t) oc 
exp[i(k n ■ x + ak 2 t)}, where k n is the wave vector of the nth beamlet and a > is a (real) 
constant. It can be checked that all the assumptions made on S are fulfilled. 

We are interested in the critical coupling X q (x, t) and its Laplacian-free counterpart 
X q (x, t) obtained by setting m~ l = in Eq. (1). These quantities are defined by 

X q (x,t) = inf{A > : (\£(x,t)\ q ) = +oo}, (5a) 
X q (x,t) = inf{A > : ( e 9*/o |S(*.t)I 2 *-) = +OQ y (5b) 

Equations (5) give the values of A at which (\S(x, t)\ g ) diverges with and without the Lapla- 
cian on the left-hand side of (1). Note that S is not assumed to be homogeneous and the 
critical coupling will depend on x in general. 



III. DISTRIBUTIONAL FORMULATION OF £(x,t;s) 



Let s be the M-dimensional Gaussian random vector the elements of which are the s n , 
and 7(2;, r) the M x M Hermitian matrix defined by 

lnm(x,r) = $*(x,r)$ m (a;,r). 

Inserting (3) into the right-hand side of (2) yields 

S(x,t;s) = [ e m^) 2 ^^M^d[x(.)l (6) 

Jx(-)eB(x,t) 

where we have made the dependence of S(x,t) on the realization of s explicit. In order 
to make (6) more appropriate to a distributional formulation it is desirable to replace the 
quadratic form s^(x(t),t)s with its monomial decomposition. One obtains 

N 



S(x,t;s)= [ e^o-M 2 ^exp X^Hs) [ Vi {x{r),T) dr 

Jx(-)€B(x,t) i=1 JO 



d[x(.)}, (7) 



with N = M 2 , and where the (fi are N real valued functions given by 7„„, \/2~Re{^ nm ), 
and \Z2Im(7 nm ), n < m. The components of the vector k(s) G M N are given by |s n | 2 , 
v / 2Re(s n s^ 1 ), and \/2\m.(s n s* m ) , n < m. It can be checked that 

ll*«ll = Ikll 2 , (8) 

where = (j2f=i h(s) 2 ^j and ||s|| = ^X^i l s i| 2 ) • We first give a heuristic 

derivation of the distributional formulation of (6). Then, we set it on a much firmer ground 
by justifying it rigorously from a mathematical point of view. 

A. Heuristics 

Inserting the identity 

N 



1 = Yi J S ~~ J ¥i( x ( T )i r ) dr^j dui, 



in the path-integral (7) and permuting the path- and M-integrals, one obtains 



S(x, t;s)= [■■■[ G x , t (u) e Xk ^ u TT dm, (9) 



with 

G x , t ( Ul ,...u N ) = [ e^^ 2dT T\s(ui - [\ z (x(r),r)dr) d[x(-)}. (10) 

Jx(-)eB(x,t) i=1 V Jo J 

As a Feynman-Kac path-integral, the expression (10) is not well defined. A possible way 
to make it meaningful consists in writing G x j as the Fourier transform w.r.t. i] of some 
function ^(x, t; rj): 

GxM = [■■■ [ t; V ) e*-" fl dm, (11) 

in which ^(x, t;rj) has a well defined meaning. Fourier transforming (10) w.r.t. u and 
permuting the path- and li-integrals, one obtains 

Jx(-)eB(x,t) 



where V(x, t; r]) is given by 



JV 



V(x, t;r]) = Y^ Vm(x, 0- ( 13 ) 



i=i 



We now observe that Eq. (12) is the path-integral solution to the Schrodinger equation 

(14) 



idt*(x, t; rj) = -±.A*(x, f ; rj) + V{x, t; rj)*(x, t; rj), 



t > 0, x E A, and 0; 77) = 1. 

This yields a well defined ^(x, t; 77). 

The permutation of path- and ordinary integrals, as well as the formal Feynman-Kac 
path-integral used in the derivation of Eqs. (9), (11), and (14) above require justification. 
The work by Cartier and DeWitt-Morette [4] suggests that we define the path-integral (6) 
by the right-hand side of (9) in which G Xjt is defined by its Fourier transform given as the 
solution to Eq. (14). We now prove the validity of this approach. 



B. The distributional formulation 

Let ty(x,t;rj) be the solution to (14) where V(x : t;rj) is given by (13) with rj E C^. Let 
at = mi x( . )eB(X:t) JjJ (pi(x(r), r) dr and k = sup a . ( . )6B(X)t) /„* (Pi(x(r), r) dr. Then the following 
lemma holds. 

Lemma 1. For every t > 0, x E A, and m E C+\{0} ; 
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(i) Gx,t defined by (11) is a distribution with compact support on M. N and suppG^j C 

[aiA] x ••• x [a N ,b N ]; 

(ii) S(x,t; s) defined by (9) is the solution to (1). 

Proof. Taking the derivative of (14) with respect to 77*, the complex conjugate of r] i: 
and using d r) *V(x,t;rj) = which follows from analyticity of V(x, t;rj) in rj [see Eq. (13)], 
one finds that d v *^(x,t;rj) evolves in time according to the same equation (14) with the 
initial condition d v *^/(x,0;r]) = 0. Thus, d v *^(x,t;r]) = for all t > and rj e C N which 
implies that *(#, t;rj) is analytic in r\. 

Let ^f(x, t; 7]) = *(:r, t; rj) exp(— it J2i=i Vi c i) where the constants q G M will be specified 
later. ^ is the solution to (14) with V given by (13) in which the are replaced by 
<Pi = ipi + Ci. Let €i = sgn[Im(r/j)]. From the equation (14) for \1> and the Schwartz inequality 
one obtains 

rf n,T,n2 _ Im(m) ||T7 ~ r||2 ^ T _ , , /• _ , = |2jd 



^11*115 = -^^HV*||2 + 2^Imfai) lji\Wd a x 



EM*) / 

i=l 

AT 

< 2||*||^ |Im(^)|sup(e^), (15) 



i=i a;eA 



4l|V*||l = -^^IIA*||? + 2Vlm(//;) / ^IWRlV.r 



i=i ■ /a 

+ 2Im^77i / [*V**] ■V^> i d d x (16) 
i=i ^ A 

AT AT 

< 2||V*||^ |Im(^)| sup(e^) + 2|| V*|| 2 ||*|| 2 J] |^| || V&IU, 



1=1 1=1 



and 



A*||2 = -^^[vA^I^Vlm^) / ^|A*|Vz 

H " ^ A 

N r 

+ 2Im^77i / A^* [*A^ + 2V*- V&l (17) 

i=i 

AT 

< 2||A^||^|Im(^)|sup(e^) 

i=i * eA 

AT 

+ 2\\AV\\ 2 J2\Vi\ (\M2\\A<Pi\\oo + 2\\V^\\V<Pi\\oo) , 



where || • || 2 and || ■ ||oo respectively denote the L 2 and uniform norms 3 on A for given t 
and r]. Both || V</3j||oo(X) and || A<^j|| OC) (t) are bounded by assumption. Integrating then the 
inequality (15) over time from to t, one obtains 



(18) 



Similarly, by integrating (16) and (17) one finds 



N 



N 



Vi 



e E,1i IM%)l/o snp xeA [ei<pi(x,T)]dT ^ 



where C\ and C<i are finite and independent of 77 and m. We now substitute (18) and (19) 
into the right-hand side of the Sobolev-type inequality below, valid for d < 3 (see e.g. [12] 
pp 106-107), 

\V(x,t;ri)\ < C 3 [||*|| 2 M) + ||A*|| 2 (i,77)] , 
with C3 finite and independent of rj and m. This yields 



\*(x,t;rj)\ < 



N / N \ : 

A + BtY,\m\ + ct 2 [Y^W 

i=l \i=l J 



3 Eili l Im (%)l So supyeAitiViiy^)} dr ? (20) 



where A, £>, and C are finite and independent of rj and to. Take 

C i = --Kl\l sup[v?i(x,r)]dr+ / inf [^(x, r)] rfr 1 , 
& I Jo xeA Jo xeA J 

and define 

Ki= [ swp[ei<pi(x,T)] dr = ^ { / sup[y?j(:r, r)] rfr - / inf [^(x, r)] dr 

JO zGA ^ I JO ^eA Jo xeA 

Note that, with this choice of q, is independent of Since k« > 0, one can bound the 
right side of (20) by 

2" 



\*(x,t;ri)\ < 



N / N \ ' 

i=l \i=l / 



(21) 



3 In the case of a vector field v(x,t) G C d , these norms arc to be understood as |H|2(t) = 
(J A v(x,t) ■ v(x,t)*d d x) 1/2 and ||u||oo(i) = sup xeA (^/v{x,t) -v(x,t)*y 
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Let G x>t be denned by Eq. (11) in which \P is replaced by By the Paley- Wiener 
theorem in the formulation given in [11] (Theorem XVI in Chapter VII), it follows from 
the analyticity of ^f(x, t;rj) in r\ and Eq. (21) that G Xtt is a distribution with compact 
support on M. N and suppG^ C [— k,i,ki] x ... x [— kn, k>n]- From the definition of \& one 
has G x j(u) = G x j{u — ct), which implies immediately by translation that G Xtt is also a 
distribution with compact support on R N and suppG^ C [«i,/3i] x ... x [(Xn,Pn\, where 
Q!j = J* * inf;,. 6 A Vi( r ) dr and = J Q * sup^gA (/9j(r) rfr. The permutation of time integral and 
space supremum (resp. infimum) can then be performed by using Lemma Al with W = ±ipi 
(see Appendix A). One obtains ctj = a« and fa — bi, yielding 

suppG^t C [ai,6i] x ... x [a N ,b N ] . (22) 

It is worth noting that, heuristically, (22) follows immediately from the formal expression 
(10) since the product of the delta functions vanishes identically outside [a±, bi] x ... x [a N , b N ]. 

It remains to prove that £(x,t;s) defined by the r.h.s. of (9) is the solution to (1). 
To this end it suffices to note that Eq. (9) can be written as £(x, t; s) = ^(x, t;r) — i\k(s)) 
which is the solution to (14) with rj = i\k(s) in the potential (13). It can be checked that the 
latter equation is indeed Equation (1) [reconstruct s^s from its monomial decomposition 
and multiply (14) by —i], which completes the proof of Lemma 1. 



IV. CONTROLLING THE GROWTH OF \£(x,t;s)\ q 

The advantage gained by recasting (6) as (9) is that the latter formulation is suitable for 
a straightforward application of the Paley- Wiener theorem (see e.g. [11] Theorem XVI in 
Chapter VII, and [6] Theorem 7.4 in Chapter VI), offering the possibility of controlling the 
growth of £(x,t;s) as ||s|| — > +oo. This is embodied in Lemma 2 below. Let s = s/\\s\\ 
be the direction of s in C M and H X;t (s) = sup x ^ eB ^ x ^ J * U(x(t), t; s) dr, with U(x, r; s) = 
52±=iH s )i ( Pi( x > T ) where k(s) = fc(s)/||A;(s)||. 

Lemma 2. For every t > 0, x G A, m G C+\{0} 7 and q a positive integer, one has 

In \£(x, t; s)\ q , TT 
limsup 1 ' n = q\H x , t (s), (23) 

||s||— »+oo ll^ll 

along every given direction s in C M . 
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Proof. From Eqs. (8) and (9), one can rewrite the left-hand side of (23) as 



,. ln\S(x,t;s)\ q 1 
hm sup j-j — j— r = q\ hm sup 



\\s — »+oo 



||fc(-)|H+oo A\K S )\\ 



In 



N 



G x , t (u)e Xk ^ u l[du, 



i=i 



Fixing the direction of s in C M also fixes the direction of k(s) in R N . Write u = vk(s) +u±, 
with u± ■ k(s) = 0, replace G Xtt (u) by its Fourier representation (11), and let r)\\ and r)± 
denote the Fourier conjugated variables of v and u±, respectively. One obtains, 

N 



[ ■ [ G X:t (u) e Afe(s) ' u TT dui = 



,\\\k(s)\\v 



(27T) 



A' 



N-l 



^t^e^i+^x)^ \\ d V±i 



i=l 



N-l 



dv Y\ du±i. 



Performing the integration over u± first, and then the one over i]±, one finds that the latter 
expression reduces to 

N 

[ .. I G x ^u)e Xk ^ u T\du t = [ g x , t (v)e x ^ v dv, 

with 



9x,t( v ) = 



(27T) 



*(x,*;77||,77 J . = 0)e iM 'iid77||, 



where ty(x,t;r)\\,r)± = 0) is the solution to (14) with V(x,t;r)) = r)\\U(x,t; s). Let a = 
inf x (.)eB(x,t) Sq U(x(t), t; s) dr and b = snp x( . )€B(x ^ | * U(x(t), t; s) dr = H Xtt (s). By Lemma 
1, g x>t is a distribution with compact support on R and suppg^t C [a, b]. This implies that 
sup{f : v G suppg Xjt } < H x j(s), and by the Paley- Wiener theorem, 

1 



limsup 

|fc(s)||-+oo A||fc(sJ 



In 



9x,t(v) e xms)llv dv 



< H Xtt (s). 



(24) 



We now prove that (24) is an equality. Suppose that 3e > such that 

1 



limsup 

||fc(s)||-+oo M\k{s) 



In 



9x,t 



v) e xms)llv dv 



< H Xtt (s) - e. 



(25) 



Then, according to the Paley- Wiener theorem, sup{t> : v G suppg^t} < H Xjt (s) — e. It is 
shown in Appendix B that sup{t> : v G suppg^} = H Xjt (s), yielding H Xjt (s) < H X}t (s) — e, 
in contradiction with e > 0. Thus, Eq. (25) is false and one obtains 

1 



limsup 



In 



9x,t 



» e x ^ s ^ v dv 



H x , t (s). 



This completes the proof of Lemma 2. 



10 



V. DETERMINATION OF X q (x,t) AND COMPARISON TO X q (x,t) 

In this section we prove the conjecture made in Ref. [2] that \ q < \ q , in the case where 
S(x,t) is given by (3). Since we wish to express the results in terms of eigenvalues of the 
correlation function (S*(x(t),t)S(x(t'),t')), we begin with a technical preliminary linking 
these eigenvalues to those of the matrix J* * 7(rr(r), r) dr. 

Let > /x 2 [#(•)] > • • • > be the eigenvalues of the covariance operator T x ^ 

acting on /(r) G L 2 (dr), defined by 

(T x(0 /)(r) = [\s*(x(r),r)S(x(r'),r'))f(r')dr', (26) 
Jo 

with < r, t' < t and x(-) G B(x ) t). Let /j(r) G l?(dr) be the eigenfunction associated 
with and define the vector cr, G C M by <7 iri = / * $*(x(r), r)/*(r) dr. From (26), (3), 

and (4) one has 

M 

(T x (.)fi){r) = y2KMr),r) ( r* m = f i t [x(-)]f l (r), (27) 



m=l 



and 



^0 

= Va jra / $;(x(r),r)$ m (x(r),r)rfr (28) 



m=l 
M r H 

lnm(x(T),T) dr 



= E 



m=l 







It follows from the last equality of (28) that any non vanishing eigenvalue of T x ^ is 
also an eigenvalue of J^(x(t),t) dr with eigenvector <7j. Conversely, any non vanish- 
ing eigenvalue of J* ^(x(t),t) dr is also eigenvalue of T x ^ with eigenfunction /j(r) = 
^ :1 Em=iC^m( :i: ( T )i r ) [ see Eq. (27)]. Thus, there is a one-to-one relationship between 
the non vanishing eigenvalues of T x (.) and jj* 7(x(r), r) dr. In the sequel /ii[x(-)] will denote 
the largest eigenvalue of T x ^ and of J* * 7(x(r), r) dr. Define 

At x ,t= sup A«i[x(-)]. 

x(-)6S(x,t) 

One can now prove the following proposition: 

Proposition. For every t > and iGA, \ q (x,t) = (q^xj)' 1 < \(x, t). 
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Proof. First we prove X q (x, t) > (q^ x ,t) 1 - Expressing U(x(t),t;s) in terms of the 
quadratic form s^j(x(t), t)s in the expression for H Xit (s), one has 



H x , t (s) = sup — 
x(-)eB(x,t) \\ s \ 



j(x(t), t) dr 



< 



Hence, by Lemma 2, 



hmsup r—r^ < q\fx X! t. 



sh+oo 



This implies that for every A < (q/J> x ,t) \ 



{\e(z,m 



d 2 s r 



e-^ 2 \£(x,t;s)\oT]^<+oo, 



(29) 



n=l 



which proves X q (x,t) > (?/4r,t) -1 - 

We now prove the inequality X q (x,t) < (qHxj)^ 1 more exactly, X q (x,t) < (q/i x j — 
+ ) -1 ]. To this end we follow the same line of reasoning as in Ref. [7]. Let A(r) = {z e 
C M : \z n \ < r, 1 < n < M}. For Eq. (29) to hold it is necessary that, for every r > 0, 

M J2„/ 



hm /••• / eHI 8+8 '"V(x,*;s + sO| 9 ft — = 0> 
'IH+ooJ J A(r) n=i n 



(30) 



along every direction s in C M . For any fixed s G C M , e~ 2s *' z,/9 £(a;, t; s + z) is an entire 
function of z E C M and hence \e~ 2s *' z / q £(x, t; s + z)\ q is subharmonic w.r.t. each component 
of z [5]. Thus, writing 



-\\s+s'\\ 2 _--||* 



2 _||„/||2 



e 11 1 e 



exp I — s* ■ s' 

q 



in the integral (30), one obtains by the subharmonicity 



/... / e-" s+s '" 2 \£(x,t;s + s')\«f]±A 
J J A{r) n=i 71 

M , 2 / 

S 'H 2 e-l s *- s 'S(x,t;s + s')\ q T] — 

IT 



A(r) 

2\~\ M -II- -' 



n=l 



> [l-exp(-r 2 )] M e-H s lH£(x,t; S )r, 



and the condition (30) implies 



lim e 



! \£(x,t;s)\* = 0, 



(31) 
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along every direction s in C M . Since every element of the matrix J* * j(x(t), t) dr is a 
continuous functional of x(-) G B(x, t) with the uniform norm on [0, t] (see the appendix B), 
its eigenvalues are also continuous functionals of x(-). Accordingly Vt > =ta e (-) G B(x,t) 
such that // Xj t — £ < fj,i[x e (-)] < /i X) t. Let <7 £ G C M (with ||<7 £ || = 1) be an eigenvector of 
f* j(x £ (t), t) dr associated with the eigenvalue ^i[x £ {-)\ and take s = a £ , then 



As) 



rt 



sup (7 £ 

x(-)6S(x,t) 
t 

7(x £ (r),r) rfr 



7(x(r),r)dr cr £ 

0"e > A*x,t - £■ 



(32) 



Thus, along the direction of a £ , Lemma 2 and Equation (32) yield 

In \S (x, t:s)\ q 



lim sup 

||s|H+oo 



and for every A > (qfx Xtt - qe) 1 , 



lim sup e 

IklH+oo 



\£(x, t; s)\ q = +oo, 



in contradiction with (31). Therefore, \ q (x,t) < (qfj> x ,t — ■ Taking e arbitrarily small 
one obtains X q (x,t) < {q^ x ,t — + ) _1 , hence X q (x,t) = (g/i Xji ) _1 . 

Finally, we always have X q (x, t) = l/q/ii[x(-) = x] (see [2]) and /i Xjt > /ii[x(-) = x] 
yields X q (x,t) < X q (x,t), which completes the proof of the proposition. 



VI. SUMMARY AND PERSPECTIVES 



In this paper, we have studied the effects of propagation on the divergence of the solution 
to a linear amplifier driven by the square of a Gaussian field. We have considered a model 
in which the propagation is that of a free Schrodinger equation with a complex mass. For 
this model, we have explicitely determined the values of the coupling constant at which the 
moments of the solution diverge. We proved that the divergence yielded by a propagation- 
free calculation, i.e. in the limit of an infinite mass, cannot occur at a smaller coupling 
constant than the one obtained with a finite mass. This extends the results of Ref. [2] where 
such an inequality was proven in the diffusion case only, i.e. imaginary mass. 

As explained in the conclusion of Ref. [2], the stumbling block to going beyond the 
purely diffusive case was to control the growth of a complex Feynman path-integral. Our 
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solution of this problem is based on the realization that, if S is given by Eq. (3), the Feynman 
path-integral can be rewritten as the Fourier integral of a distribution with compact support 
(Lemma 1). Control can then be obtained as a consequence of the Paley- Wiener theorem 
(Lemma 2). 

In conclusion we outline some possible generalizations of this work. From a practical 
point of view, it would be interesting to find out whether there exists a class of S of the 
form (3) for which there is no propagation effects on the onset of the divergence, i.e. for 
which X q (x,t) = X q (x,t). In addition, since most Gaussian fields of physical interest admit 
Karhunen-Loeve-type expansions, it would also be very interesting to find a way to generalize 
our solution to the case where the finite sum (3) is replaced by an infinite sum. Other prob- 
lems involve relaxing some of the assumptions in (1). For instance, under what conditions 
on S do the results carry over to the case where A is replaced by M. d and £(x, 0) G L 2 (R d ). 
It should also be checked whether our solution of the problem is robust with respect to the 
initial condition. If the answer is no, the size of the set of £(x, 0) for which our results do not 
hold should be estimated according to physically relevant measures on the space of £(x, 0). 
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APPENDIX A: ON THE PERMUTATION OF TIME INTEGRAL AND SPACE 
SUPREMUM 

This appendix is devoted to the proof of the following lemma, 

Lemma Al. Let A be a compact pathwise connected metric space [with distance de- 
noted by d(-,-)J, and t > a real number. Let W be a real function on A x [0,t] such 
that W(-,t) is continuous in x uniformly in r G [0,t], and Vrr G A, W(x,-) is piecewise 
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continuous with a finite number of discontinuities. Then, for any 16A 

/ sup W(y, r) dr = sup / W(x(t),t) dT, 
Jo yeA x(-)eB(x,t) Jo 

where B(x,t) is the set of continuous paths in A satisfying x(t) = x. 

Proof. We obviously have 

/ sup W(y, t) dr > sup / W(x(t), t) dr, 
Jo yeA x(-)eB(x,t) Jo 

and it remains to prove the inequality in the other direction. First we assume W G C°(A x 
[0, t}). Since Ax [0, t] is compact, W is uniformly continuous, and for every e > we can find a 
number 5 = 5(e) > such that if max{d(x, x'), \t—t'\} < 5(e), then \W(x, t)—W(x', t')\ < e. 
Moreover, 5(e) tends to zero with e. If t/5(e) is not an integer, it is convenient to replace 
5(e) by the smaller quantity t/(l + [t/5(e)]) where [•] denotes the integer part, and from now 
on we will assume that t/5(e) is an integer. For a fixed e > 0, let = N(e) = t/5(e), and 
let 1Z e be a finite partition of A by sets of diameter at most 5(e). We observe that for any 
< q < N - 1 one has 

SUp SUp OsC Fx [ qt /N,(q+l)t/N]W < e, (Al) 
F&TZ e 0<q<N-l 

where "Osc" denotes the oscillation of the function (namely, its sup minus its inf). We now 
choose once and for all a point (xp, q , T q ) in each F x [qt/N, (q + l)t/N]. 
For each < q < N - 1 and any r G [qt/N, (q + l)t/N] one has 

supW(y,r)< sup sup W 

yeA FeTZ e Fx[qt/N,(q+l)t/N] 



< sup 

F£K C 



sup W + W(x Fq ,T q ) - inf W 

Fx[qt/N,(q+l)t/N] ' Fx[qt/N, (q+l)t/N] 



< sup W(x F>q ,T q ) + sup OsC Fx[qt / N) ( q+ i) t /N]W 

F£Ke F£Ke 

< SUp W(x F:q ,T q )+ SUp SUp OsC Fx [qt/N , ( q+ l)t/N] W 
FeTZe FGllz 0<q<N-l 

< sup W(x F>q , T q ) + e, 

where we have used the inequality (Al). Thus, choosing an atom F q G 1Z e such that 
su Pfg7^ W(x F:q , T q ) = W(x Fq , q , T q ), one has for any r G [qt/N, (q + l)t/N] 

supW(y,T)<W(x Fq , q ,T q ) + e. (A2) 
yeA 
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Let yi and y 2 be two given points in A. We now define a continuous path x e in A from y\ to 
y 2 . For any 1 < j < N — 1 we choose a family of continuous paths Xj from [— t/2N 2 , t/2N 2 ] 
to A satisfying Xj(—t/2N 2 ) = xp^j-i and Xj(t/2N 2 ) = x Fj: j. We also choose a continuous 
path x from [0,t/2N 2 ] to A such that x (0) = y\ and x (t/2N 2 ) = xf 0j0 , and a continuous 
path x N from [—t/2N 2 ,0] to A such that x N (0) = y 2 and x N (—t/2N 2 ) = xf n _ u n-i- The 
continuous path x e is defined by 

' x (t) for < t < t/2N 2 , 

x{t) = < x F q , q for qt/N + t/2N 2 < r < (q + l)t/iV - t/2N 2 , 

I x 9 (r - gt/7V) for qt/N - t/2N 2 < r < qt/N + t/2N 2 for q ^ 0, 

k x n (t - t) for * - t/2N 2 <r<t. 

One can observe that the Lebesgue measure of the time domain over which x e (r) 7^ xp q ^ q 
is at most equal to t/N. Since A is compact and Vr G [0,t], W(-,t) G C°(A), there is a 
finite number M > such that for any (x,r) G A x [0,t] one has \W{x, r)| < M. Thus, 
\W(x Fq , q , r) - W(a; e (r), r)| < 2M for any r in [qt/N, (q + l)t/N\. Using the latter estimate, 
the remark below (A3), and (A2) one obtains 

rt N-l ,( q+ l)t/N 

I sup W(y, t) dr = \^ / SU P W{y, r) dr 

Jo 2/eA J qt/N yeA 

N-l Mq+l)t/N 

<et + J2 W(x Fqtq ,r)dr (A4) 

q=0 J H/ N 

2Mt f f 
< et + — + J o W(x e (r),r)dr, 

Now, assume that there is a finite set of times independent of x, {t±, tl}, such that 
Vrr G A the set of times at which W(x, •) is discontinuous is a subset of {t±, ...,tl}- Thus, 
(A4) applies in each time interval [r^Tj+i], < i < L, with r = and r^+i = t. Let 
y , yi, yL+i be L + 2 points in A with yi+i = x. Let x t be a continuous path in A passing 
by x e (Ti) = yi and defined by (A3) in each time interval [Tj,r i+1 ]. From Equation (A4) in 
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which one writes Cj(e) = e + 2M/iVj(e), it follows 





i(-)eB(i,t) 7o 



where the last inequality results from the fact that x e G B(x,t). The proof of Lemma Al 
for the class of W considered in this paragraph is completed by taking the limit e — > and 
observing that the Cj(e) tend to zero with e. 

We are now ready to prove Lemma 1 in the general case. Since W(-,t) is continuous 
in x uniformly in r G [0, t], for any e > 0, we can find 5 = 5(e) > such that if d(x, x') < 5, 
then 



Since A is compact, we can find a finite covering by open balls of radius at most 5/2, and 
therefore a finite partition of unity (xk) by continuous functions whose support has diameter 
at most 5 (see [3], paragraph 4.3 in Chapter IX). For any k we choose once for all a point 
Xk G suppxfc, and define the function 



For each fixed r, this function is obviously continuous in x, and for fixed x, it is piecewise 
continuous in r, with the possible discontinuity points belonging to a finite set which can 
be chosen independent of x. From the previous result we have 



sup \W(x,t)-W(x',t)\ < e. 
re[o,t] 



W e (x,r)=Y J W(x k ,r) Xk (x). 



k 






sup W(y, r) dr < et + I sup W e (y, r) dr 




= et + sup / 

x(-)eB(x,t) Jo 




< 2et+ sup / W(x(: 

x(-)eB(x,t) Jo 

Since the estimate holds for any e > the general result follows. 
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APPENDIX B: DETERMINATION OF THE SUPPORT OF 

(x,l ) 

Let g£f(u) be a distribution with compact support on K. whose Fourier transform, 
^ (m \x,t;r]) = {Fg { ™t){r]) wi th G M, is the solution to (14) with V(x,t;rj) = rjU(x,t;s), 
where U(x,t; s) = k(s)i(fi(x,t). This appendix is devoted to the determination of the 
support of g^ly We have modified the notation used in the text to make the dependence 
on m explicit. 

We begin with a technical lemma that will be useful in the sequel. Let C^°(R) denote 
the set of all smooth compactly supported functions in R, 

Lemma Bl. For every t > 0, x G A, and f E Cq°(M.), f R g^l\v)f(v)dv is an 
analytic function of m on C + = {m E C : Im(m) > 0} ; and f R g^\v)f(v)dv = 
lim 7 ^ + J R g^} +%1 \ v ) f ( v ) dv for each real m/0. 

Proof. As a Fourier transform of a function with compact supports on R, (J 7 f)(rj) is 
an analytic function of 77 e C. We have seen at the beginning of the proof of Lemma 1 
that ty( m ^(x,t;r]), with m G C + , is also analytic in rj. Furthermore, if A is a torus and V 
is bounded on A (which is the case), then (i) \&( m )(x, t\rj) is analytic in m E C + ; and (ii) 
Vr] E C, ^ m \x,t;rf) = lim 7 ^ + \I/( m+i7 )(a;, t; rj) for each real m/0. We are indebted to K. 
Yajima for the proof of the latter result that we reproduce here for the sake of completeness 
[13]. Define 

U m (t) = exp (i^J , t > 0, m E C + , 
and write the initial value problem (14) in the the form of integral equation 

#M(t) = l-i v [ U m {t - T)V{T)¥ m \T) dr. 
Jo 

Here V(t) is the multiplication operator with U(x,t; s). Let B denote the space of bounded 
operators in L 2 (A). By Fourier series expansion, it is evident that (a) \\U m (t)^ m \t)\\ 2 < 
||^( m )(t)|| 2 , viz. U m {t) E B and \\U m (t)\\ < 1; (b) the function [0, 00) x C+ 3 (t,m) -> 
U m (t) E B is strongly continuous [viz. (t,m) — > limit) f E L 2 (A) is continuous for every 
/ G L 2 (A)]; and (c) for every t > 0, m — > W m (i) G i3 is analytic for m G C + and {d/dm)U m {t) 
is norm continuous w.r.t. (t, m) G [0, 00) x C + . It follows from the boundedness of V that 



18 



the Dyson expansion [8] 

D m (t) = U m (t) -it] f U m (t - r)V{r)U m {r) dr + ■ ■ ■ + 
Jo 

(-ivT [ U m {t - r n )V(r n ) ■ ■ ■ Vir^Umin) dr ± ■ ■ ■ dr n + ■■■ 

J0<Tl<-<T n <t 

converges in the operator norm of B uniformly w.r.t. (t,m) in every compact subset of 
[0,oo) x (C + \{0}). Thus, the operator D m (t) enjoys the same properties (b) and (c) men- 
tioned above as an operator valued function of t and m. It is easy to check that D m (t) defines 
the propagator for (14) and is unitary if m is real. Hence the solution to (14) satisfies the 
properties (i) and (ii). 

Analyticity of ^ (m) (x, t; rj) and (J 7 /) (77) in m E C + and 77 E C implies that 



g^\u)f(u)du= [ ¥ m \x,t- v ){Ff)(- v ) 

Jr 



drj 
2n 



is an analytic function of m E C + . By Eq. (20) and the fact that none of the constants 
A, B, and C depend on m, \^( m \x,t\rj)(J- f)(—rj)\ is bounded by an integrable function 
of rj independent of m, from which it follows that the m-limit and the ?7-integral can be 
interchanged. Thus, according to (ii), one finds that for each real m/0 



lim + / g { Z +il \u)f{u)du = lim / n)(rf)(-rj)p- 

7^0+ J R 7 -*0+ J R ATT 



= [ \hnjF^\x,t; V ){Ff)(-Ti)p- 

Jr 7^0+ 2tT 

= [ ¥ m \x,t;r ] )^f)(-r ] )^- = [ g^(u)f(u)du, 

which completes the proof of Lemma Bl. 

Let a = mi x{ . )eB{Xjt) J*U(x(r),r;s)dr and b = sup x{ . )eB{X:t) J* U(x(r), r; s) dr = 
H x ,t(s)- One has the following Lemma: 

Lemma B2. For every t > 0, x E A, and m E C+\{0} ; the support of g x ^ is equal to 

[a,b\. 



Proof. First, consider the case m = ry, 7 > 0. Denote by a[x(-)] the functional 
•)] = f*U(x(T),r;s) 
Vx, y E A and < r < t, 



a[x(-)] = J* U(x(t),t; s) dr. By boundedness of W over A x [0,t], 3A > such that 



\U(x,t;s) -U(y,r;s) 



f VU(x', r; s) ■ dx' 



< A\\x-y\\, 
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where x, y denotes the segment of geodesic from x to y and || ■ || is the usual Euclidean 
distance. This implies that Vx(-),y(-) G B(x,t), 

\a[x(-)]-a[y(-)]\ < f \U(x(r),r; s) - U(y(r),r; s)\ dr 
Jo 

< A I \\x{r)-y{r)\\dT 

Jo 

< At sup \\x(t) - y(r)||, 

0<T<t 

which shows that ct[x(-)] is a continuous functional of x(-) G B(x, t) with the uniform norm on 
[0, t}. Let h G C£°(1R) a real positive test function with support in [a, b] and sup uGR h(u) = 1. 
From the continuity of it follows that Ete (-) G B(x,t) such that /i(a[xo(-)]) = 1. By 

continuity of h and h follows that Ve > 0, 35 > such that |/i(o;[x(-)]) — 1| < e for 

every x(-) G B (5) = {x(-) G B(x,t) : sup 0<T<t ||x(r) — ^o( r )|| < £}• Take e = 1/2, in this 
case h(a[x(-)]) > 1/2 for every x(-) G S (5) and one has 

/ g%\u)h(u)du = [ e-^o^) 2 ^M«K-)])^(-)] 

</R Jx(-)eB(x,t) 

> [ e-*ti ±{T)2dT h(a[x(-)\)d[x{-)\ 

Jx( : )£B (8) 

e -i/o*W 2 *d[a;(.)]. 

--(■)GBo(5) 

Since the set of the Brownian paths #(•) that are in B (S) has a strictly positive Wiener 
measure, the last term is strictly positive and one finds 

/ g ( :j\u)h(u)du>0. (Bl) 
Jul 

If there was an open subset of [a, b] not intersecting the support of g^j\ it would be possible to 
choose the support of h outside the one of g%j\ yielding J g^j\u)h(u) du = in contradiction 
with Eq. (Bl). Thus, for every x G A and 7 > 0, the support of g^J* is equal to [a, b]. 

Consider now the general case m G C+\{0}. Assume that there is an open subset of 
[a, b] not intersecting the support of g^f ■ In this case it is possible to choose the support of 
h outside the one of , yielding J g^\u)h(u) du = 0. By Lemma Bl the support of g^} 
must vary continuously with z G C+, whence the support of h can be taken small enough 
such that there is a open subset V(m) C C + with m G V(m) and J g^J(u)h(u) du = 
identically in V(m). From the analyticity of J g^J{u)h(u) du in z on C + (Lemma Bl), it 
follows immediately that J g^J(u)h(u) du = identically in all C + , in contradiction with 
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Eq. (Bl), which completes the proof of Lemma B2. In particular, one has sup{t> : v G 
suppy^} = b = H Xtt (s), which is the result used in the proof of Lemma 2. 
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